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We derive all second post-Newtonian (2PN), non-precessional effects of spin-orbit coupling on 
the gravitational wave forms emitted by an inspiraling binary composed of spinning, compact bodies 
' in a quasicircular orbit. Previous post-Newtonian calculations of spin-orbit effects (at 1.5PN order) 

, relied on a fluid description of the spinning bodies. We simplify the calculations by introducing into 

post-Newtonian theory a 5-function description of the influence of the spins on the bodies' energy- 
momentum tensor. This description was recently used by Mino, Shibata, and Tanaka (MST) in 
Teukolsky-formalism analyses of particles orbiting massive black holes, and is based on prior work 
by Dixon. We compute the 2PN contributions to the wave forms by combining the MST energy- 
momentum tensor with the formalism of Blanchet, Damour, and Iyer for evaluating the binary's 
radiative multipoles, and with the well-known 1.5PN order equations of motion for the binary. Our 
results contribute at 2PN order only to the amplitudes of the wave forms. The secular evolution of 
the wave forms' phase — the quantity most accurately measurable by LIGO — is not affected by our 
. results until 2.5PN order, at which point other spin-orbit effects also come into play. We plan to 

evaluate the entire 2.5PN spin-orbit contribution to the secular phase evolution in a future paper, 
CO . using the techniques of this paper. 



*q ' I. INTRODUCTION 

cr 

$H . Inspiraling compact binaries are one of the main classes of gravitational wave source to be targeted by the com- 
ing generation of ground-based laser interferometers such as LIGO, VIRGO, GEO600, and TAMA (lj. There are 
, two reasons for this. First, binary coalescences are expected to occur fairly often within detection range of "en- 
■ hanced" interferometers Astronomical lore estimates several neutron-star-neutron-star coalescences per year 
within 200 Mpc HU and a similar rate of black-hole-black-hole coalescences within 200 Mpc to 1 Gpc . Second, 
the signal from the final moments of inspiral is characterized by a complicated phase evolution containing detailed 
information about the physical parameters of the binary, such as the masses of the bodies and their spins about their 
own axes |TJ. 

Because inspiral signals have such a complicated structure, and because they last many cycles within the frequency 
bands of ground-based interferometers, they are ideal candidates for the use of matched filtering [Q . Matched filtering, 
a signal-processing technique well-studied in the context of radar, can be used both to search for signals in noisy 
data and to estimate parameters once a signal is found. Matched filtering essentially entails cross-correlating noisy 
interferometer data with a set of theoretical template wave forms. If a template wave form is a good approximation 
to the signal wave form, the cross-correlation enhances the signal-to-noise ratio. In the context of matched filtering, a 
good approximation means (roughly speaking) one in which the phase evolution of the template matches that of the 
signal to within a half cycle out of the total spent in an interferometer's band. Because signals are expected to last 
up to tens of thousands of cycles in the bands of some interferometers, the templates must match any possible signal 
to a correspondingly high degree of precision. 
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Currently there is no exact solution to the generic two-body problem in general relativity. Thus, inspiral wave-form 
templates are constructed using approximation schemes which must be carried out to high precision to be useful for 
matched filtering. These approximation schemes can be broadly grouped into two categories, the post-Newtonian 
approach and the black-hole perturbation approach. 

The post-Newtonian approach is the longtime standard for gravitational wave generation. It involves expanding the 
Einstein equations and equations of motion in powers of the binary's orbital velocity v/c and gravitational potential 
GM/rc 2 ~ (v/c) 2 , where the order in GM/rc 2 is referred to as the post-Newtonian (PN) order. Concurrently, the 
gravitational wave forms and luminosity are expanded in terms of time derivatives of symmetric, trace-free (STF) 
radiative multipoles, which are expressed as integrals of the matter source and gravitational fields. The radiative 
multipoles are combined with post-Newtonian equations of motion to yield explicit expressions for the wave forms, 
including the secular evolution of orbital phase and frequency due to radiation reaction. Recent summaries of the 
two main versions of the post-Newtonian approach are given by Blanchet, Damour, and Iyer || and by Will and 
Wiseman 0. The post-Newtonian expansion of the wave forms of a binary currently has been carried out to 2.5PN 
order jlO| [i.e., to 0(v/c) 5 beyond leading order radiation reaction and 0(v/c) 5 beyond Newtonian gravity in all 
non-radiation reaction effects] in the case where the two bodies orbit in a quasicircular fashion and do not spin about 
their own axes. (By "quasicircular" we mean orbits that are circular aside from gradual inspiral due to gravitational 
radiation reaction.) 

In the case where the bodies do spin, there are three types of spin effects to be considered. Effects of the first 
type, due to precession of the plane of the orbit, modulate the amplitude and frequency of the gravitational radiation 
in a complicated, non-monotonic fashion. Secular or dissipative effects, due to radiation reaction, contribute to 
the (monotonic) phase and frequency evolution of the orbit; and non-dissipative effects contribute directly to the 
amplitudes of the various harmonics of orbital frequency in the wave forms, without affecting their phase evolution. 
All three types of effects can be further divided into spin-orbit contributions (i.e., terms involving one spin only) 
and spin-spin contributions (interactions between spins). Precessional effects were first extensively investigated by 
Apostolatos et al. and by Kidder jl2| and were found to complicate matters considerably. Therefore, like most 
other treatments of spin, ours will investigate the case where there is no precession — i.e., the spins are parallel or 
antiparallel to the orbital angular momentum — leaving precession for future studies. 

Non-precessional spin effects have been evaluated by Kidder, Will, and Wiseman 13 lj] only to lowest order: 1.5PN 
for dissipative (and 1PN for non-dissipative) spin-orbit effects, and 2PN for spin-spin effects.^ The main reason for 
the discrepancy in progress between the spinning and nonspinning cases is the form of the matter source used in the 
Einstein equations. In the nonspinning case it is simple to write the energy-momentum tensor as a Dirac 5-function, 
which greatly simplifies the calculations. In order to derive spin effects, Refs. Jl^Jl^] treated the bodies as uniformly 
rotating balls of perfect fluid. The perfect fluid energy-momentum tensor was integrated over a finite spatial volume, 
which made the multipole integrals much more cumbersome than in the 5-function case and introduced additional 
complications in the definition of the binary's center of mass. The net result was that spin calculations at a given 
post-Newtonian order seemed to require as much effort as spinless calculations at higher post-Newtonian order, and 
spin calculations were not pursued any further with this approach. 

The more recent black-hole perturbation approach obtains high-order (in some cases exact) expressions for the 
influence of radiation reaction on the orbital phase which are valid in the limit of extreme mass ratio. The basis of 
this approach is the perturbation of known, exact solutions of the Einstein equations (the Schwarzschild and Kerr 
spacetimes) with a test body using the Teukolsky equation [l4| or an equivalent. During the last several years, 
analytical techniques for post-Newtonian expansion in the context of the black-hole perturbation approach have been 
developed to very high orders in v/c (for recent reviews, see fl5|l). However, most black- hole perturbation papers 
treat the test body as a nonspinning point particle with a 5-function energy-momentum tensor, and thus do not give 
results for the case of two spinning bodies. 

Recently, the black-hole perturbation approach has been extended to the case of two spinning bodies 
In pH) , Mino, Shibata, and Tanaka calculated the gravitational wave forms and radiation reaction of a spinning 
particle falling into a Kerr black hole. In Tanaka et al. obtained an expression for the non-precessional 2.5PN 
spin-orbit contribution to the secular phase evolution of a binary composed of a spinning test particle in quasicircular 
orbit around a Kerr black hole. These results were obtained using an energy-momentum tensor for the test body 
which mimics the effects of an extended, spinning object but can be expressed in terms of a (5-function for ease of 



1 Like non-spin effects, spin effects appear in the secular phase evolution of the wave forms at a certain order, and at every 
order in v/c (0.5PN order) beyond it except for the first. Thus, spin-orbit effects appear at 1.5PN, 2.5PN, 3PN. . . orders and 
spin-spin effects appear at 2PN, 3PN, 3.5PN. . . orders. 
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calculation. This "spinning-particle" (5-function ener gy- m omentum tensor is based on the work of Dixon |18|. We call 
it the MST tensor after Mino, Shibata, and Tanaka jL6|, who distilled it into the compact form we will use. 

In this paper, we use the MST energy-momentum tensor for the first time in the curved-space, post-Newtonian 
approach to derive new gravitational- wave generation results. (Cho fl9[| , in work parallel to our own, has recently used 
a similar approach to re-derive the wave forms of Kidder, Will, and Wiseman (l3) in a slightly different form.) We 
reproduce (with a shorter calculation) the 1PN and 1.5PN spin-orbit corrections to the radiative multipoles derived 
in We also derive all of the (previously unknown) 2PN non-precessional spin-orbit corrections to the wave forms, 
by calculating 2PN spin-orbit corrections to the radiative multipoles and combining them with the well-known 1.5PN 
equations of motion (in which there is no 2PN spin-orbit term). Because of the harmonics of the orbital frequency 
involved, there is no 2PN spin-orbit contribution to the radiation reaction-induced secular phase evolution of the wave 
forms (the most accurately measurable effect). 

In the future, we plan to use the methods of this paper to calculate all the non-precessional 2.5PN spin-orbit effects, 
including the nonvanishing radiation reaction and resulting secular evolution of the frequency and phase of the wave 
forms. That secular evolution is likely to be quite important for data analysis. Investigations by Tagoshi et al. p0[ , 
comparing post-Newtonian expansions to exact numerical results in the test-mass limit, indicate that spin effects are 
important for extraction of information from observed waves at least up through 3PN order. 

This paper is or gan ized as follows. In Sec. |n] we present the MST energy-momentum tensor [^6| and review its 
properties. In Sec. Ill we review the post-Newtonian expansions of basic variables used in our calculations. Then in 
Sec. |y| we calculate the STF radiative multipoles needed to obtain the 2PN spin-orbit terms in the wave forms. In 
Sec. M we evaluate all the 2PN (non-precessional) spin-orbit terms in the wave forms of a binary in quasicircular orbit 



with spins parallel or antiparallel to the orbital angular momentum, and in Sec. VI we briefly discuss their significance. 
In an Appendix we use our methods to derive the 1PN and 1.5PN STF radiative multipoles, and compare with the 
results of Refs. [||Jl|. 

Throughout this paper, we use units such that Newton's gravitational constant and the speed of light equal unity. 
We also use the tensor notation conventions of ||B|]: curved brackets () on tensor indices to indicate symmetrization, 
square brackets [] to indicate antisymmetrization, and angled brackets () or the superscript STF to indicate the 
symmetric trace-free part. A capitalized superscript L indicates a multi-index i\ ■ ■ ■ ig; e.g., I L represents J* J in the 



case £ = 3. We also write outer products of vectors in shorthand; e.g. 
indices run from to 3, and Latin indices from 1 to 3. 



-ijk 



and 



Greek 



II. SPINNING PARTICLE ENERGY-MOMENTUM TENSOR 

Our starting point is the spinning particle energy-momentum tensor given in terms of the Dirac (5-function ]16| 



T al3 (x) 



0)i 



(x-z(t)) 



S^ a (x,T)u P \x,T) 



<t)) 



(2.1) 



Here z' 1 (t) is the worldlinc of the particle, u M (r) = dz^/dr, p^{t) is the particle's linear momentum, and S^ v {t) is an 
antisymmetric tensor representing the particle's spin angular momentum. We focus only on spin-orbit interactions, 
i.e. discard all terms higher than first order in spin. In this case r becomes the particle's proper time and u* 1 becomes 
its four- velocity [see Eq. (2.4) of Ref. fu|]. 

The bitensors p a (x, t), u a (x, r), and S a ^(x, r) are spacetime extensions of u^, and S^ v away from the particle's 
worldlinej^j defined by 



u a (x,r)^g^(x,z(r))u^T), 
S^(x, t) = g^{x, z(r))gi(x, z(r))S^(r). 

Here g a M (x, z) is a bitensor of parallel displacement with the properties 



(2.2a) 
(2.2b) 
(2.2c) 



2 We use indices a, /3, 
the worldline z(t). 



to denote quantities associated with the field point x, and /i, 



to denote those associated with 
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lim gl(x,z(T))=6Z, 
lim Vpg a li (x,z(r))=0. 



(2.3a) 
(2.3b) 



The definition of is arbitrary up to the choice of a spin supplementary condition (the analog of a gauge 
condition). We use 



S^X = 0. 



(2.4) 



Note that in post-Newtonian theory at least three spin supplementary conditions are in common use. We choose (2.4) 
because it makes our radiative multipoles consistent with the standard post-Newtonian equations of motion, thus 
simplifying the calculations (cf. Ref. |l2[ , Appendix A). We introduce a spin vector S 7 ' which is related to the spin 
tensor by 



= 0, 



(2.5a) 
(2.5b) 



where e^ p<T is the Levi-Civita tensor.^ The s pin s upplementary condition is identically satisfied by ( |2.5a ). On the 
other hand, we need to impose the condition (2.5t) on to fix the one remaining degree of freedom S". 

For later convenience, we separate the MST energy-momentum tensor fl2.l| ) into the usual point-particle piece (the 
first term) plus a spin-orbit piece, 



8^(x-z(t)) 



-(so)(- T ) = - / drV 7 

When evaluating the radiative multipoles of a system of masses, we encounter integrals of the form 

d 3 xF L (x)T a(3 (x). 



(2.6) 



(2.7) 



The spin-orbit contribution to this expression can be evaluated by substituting (|2.6| ), using Leibniz' rule to rewrite 
the integral, and discarding the spatial integral of a three-divergence. The result for a many-body system is 



f d 3 xF L (x)T«%(x)=Y, 



7 (a 0) ^7 



d^F L 



- Oo ( ^ 



7 C" 5 !",' 3 ) I c7(a F /3) „,5 



1 1 8 V A 



F 



L 1 



(2.8) 



where A labels the bodies, v a = u a /u°, and <9 7 is shorthand for d/dx 1 evaluated at x = xa- 



III. POST-NEWTONIAN EXPANSIONS OF BASIC VARIABLES 



We now switch from fully covariant expressions to post-Newtonian expansions in harmonic coordinates. Spatial 
indices on the right hand sides of the equations in this section can be raised and lowered freely with the Kronecker S. 



We use the expansion parameter e which is related to the orbital variables by e ~ M/r 



where M is the total 



mass of the system, r is the orbital separation, and v the orbital velocity. We assume the bodies are compact, i.e. 
each body's spin has magnitude \Sa\ ~ X m A' wnere X is °f order unity (see Q for further discussion). 

When evaluating the post-Newtonian expansions of basic variables in this section and the radiative multipoles in 
Sec. IV, we encounter divergent expressions — in our case, self-interaction terms. Such divergences arc inevitable when 
using any (^-function source, and we follow previous authors in discarding them (see the discussion at the end of 
Ref. jq], Sec. II). We do not claim any rigorous justification for doing so; however, since it is asserted in the non- 
spinning case H that this procedure can be justified to 0(e 2 ), and since we consider corrections only up to 0(e) 
beyond lowest order spin effects, we expect that the formal use of the (^-function is justified to the same degree as in 
the non-spinning case. Informally, we note that the usual post-Newtonian equations of motion for spinning bodies 



We use the convention denned in Eq. (8.10) of Ref. 121 
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can be obtained by taking the divergence of the MST energy-momentum tensor (2.1) and discarding self-interaction 



divergences 16 



The metric components in harmonic coordinates are well known j9| as 

. 90 = -[1-2C/ + O(e 2 )], 

g tJ =5 tJ [l + 2U + 0{e 2 )}, 
sf=g = l + 2[/ + 0(e 3 / 2 ), 

where only the lowest-order expression for the potential U is needed, 

m A 



u(x) = ]T 



X - X A \ 



+ 0(e 2 ). 



(3.1a) 
(3.1b) 
(3.1c) 
(3.1d) 



(3.2) 



By differentiating ( p.l|) , we find the dominant Christoffel symbols 

1 io — 1 oo — " ! 

T) k = 8 ij a k + 5 ik a j -8 jk a\ 



(3.3a) 
(3.3b) 



where a 1 — dJJ . All oth ers a re of higher post-Newtonian order, and can be neglected for the purposes of this paper. 
The metric components (3.1), together with the condition u^u^ = — 1, give us the expansion of the four- velocity 



!+ — + [/ +0(e 2 ) 



u = v 



1 + I y + u ) + °( e 



(3.4a) 
(3.4b) 



We e xpres s t he c ompo nent s of the spin tensor in terms of the spatial components of the spin vector by combining 
(f2~5a]), j2~5b|), (IO), and fO) to obtain 



;>i0 



(« x Sf + 0{e 3/2 ) 



+0(e 2 ), 



(1 + _ w 2 _ _ . 



(3.5a) 



(3.5b) 



where e ijfc is from here on used to indicate the antisymmetric symbol [ijk]. Substituting ( |3.5| ) back into (2.6), we find 
that the post-Newtonian orders of the components of ^Vgo) are 



TnOO rpij 

1 (SO) ~ J (SO) 



,0 



mv 3/2\ 

m ~ m x 0(e J/z ), 

r 
m 



m — ~ m x 0(e), 
r 



' (SO) 

where m is the mass of either body. This contrasts with the point-mass order counting, 

i (PM) ~ TO ' 

T»° M) ~ m x 0(e 1 /2 ); 
T^ M) ~ m x 0(e). 



(3.6a) 
(3.6b) 

(3.7a) 
(3.7b) 
(3.7c) 



We also note that the second and third terms on the right hand side of ( |2.8| ) are O(e) with respect to the first if F L 
is an outer product of position vectors (as is the case when computing multipoles). 
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IV. STF RADIATIVE MULTIPOLES 



In this section we calculate the symmetric, tr ace-fr ee (STF) radiative multipolcs necessary to obtain the 2PN 
spin-orbit contributions to the wave form. In Sec. IV C| we specialize to the case of non-precessing orbits. 

The STF radiative multipoles are given to 0(e 5 ' 2 ) by Blanchet [Eq. (4.3) of Ref. @]. However, we only need the 
O(e) expressions 



I L (t) 



J L (t) 



d 3 xx L (T 00 +T") 
1 d 2 



4(2^ + 1) d 



d 3 xx La T 0a 



2(2£ + 3) dt 2 



ab{i 



d 3 xx L ' 1)a {l 

I d 2 



(£ + l)(2i + 3) dt 
d 3 xx L \x\ 2 (T 00 + T u ) + 0(e 2 ), 
21 + 1 d 



■Ob 



2(2£ + 3) dt 2 



d x x 



L-l)a 



xl 2 T 0b 



{£ + 2){2£ + 3) dt 
0(e% 



d 3 xx L ~ 1)ac T bc 



where x L denotes the symmetric trace- free part of x L . We also define I^o) ano ~ ^{SO) substituting T ( 



(SO) 



(4.1) 



(4.2) 



for T» v in 



(4.1) and (4.2). In ( f4.l|) and (4^2) we have discarded self-interaction terms, which are always divergent when using a 
(5-function source. Wc have also discarded terms involving only gravitational potentials (referred to as "non-compact" 
terms in ||) , whose spin-orbit contributions do not appear until higher post-Newtonian orders than considered in this 
paper. 

Spin-orbit corrections to the multipoles and wave form follow an order-counting scheme different from the usual 
point-mass terms. Substituting ( |3.6| ) into ( |4.1| ) an d (f4.2| ), we find that the lowest order spin-orbit correction to a 
multipole appears in the current quadrupole |12| . This term contributes to the wave-form amplitude at 1PN 

order, but because it is a subharmonic of the dominant (mass quadrupole) radiation, it does not contribute to the 
radiation reaction until 1.5PN order. The next-order effects appear in the current octupole J^g ) an d the mass 
quadrupole i?g ), an( i contribute to the wave forms and radiation reaction at 1.5PN order. (These are the terms 
given in we evaluate them with our methods in the appendix.) Following this progression, the 2PN wave forms 
require evaluation of J/gQ) an d ^(so) ^° l° wes t order, and of J%q-\ to 0(e) beyond lowest order. 



A. TV-body case 



We first evaluate the spin-orbit contributions to the multipole integrals (4.1) and (4.2) as sums over N bodies. 
The expression for the curre nt h exadecapole is the easiest to evaluate. It is needed only to lowest order and thus 
involves only the first term in (4.2), 



ijk£ 
(SO) 



e Qb < 1 / d 3 x X 3 k ^ a T^ 0) , 



which is straightforwardly obtained from the first term of (2.- 

5 

A 

The mass octupole is also needed only to lowest order, 



jijkt _ 5 ( Oi „jki\ 

J (SO) - 2 \°A X A ) 



STF 



J-ijk _ / j3 ijk ( rpOQ , rpaa 
1 (SO) — a xx \/(SO) + J (SO) 



7_d_ 

9dt 



d 3 xx l: > ka T^ 0) . 



(4.3) 



(4.4) 



(4.5) 



Again, using the first term of (2.8) it is straightforward to evaluate the integrals. When evaluating the time derivative 
in the second term we neglect time derivatives of the spins (i.e., precession). We could do this even if considering spin 
precession, because those derivatives appear 0(e) beyond the spins themselves [cf. Eqs. (F18,F19) of Ref. g, where 



G 



due to a typographical error a factor of (M/r) 3 was omitted from in front of the last term in each equation] 
left with 



.Q We 



are 



jijk _ \ 
1 (SO) ~ 



n STF 



"-{v A x SaY^a - x S A ) l x A v A 



(4.6) 



Because the two-index current moment is needed to 0(e) we must keep all three terms in (4.2), 



jij _ ab{i 
"'(SO) — fc 



1 J- 

j)acrpbc , 1 u I ,3 j)arpOb 



(SO) 



Again, time derivatives of the spins appear only in higher-order terms and may be discarded at this order, 
evaluating the integrals according to (2.8), we obtain 



(4.7) 
Carefully 



TV 

J {SO) 



E 



i qj 
2 X A°A 



1 

14 



VA ■ X A V l A S A + -Sa ■ X A v 'l + 2gXA ■ XAd A S A — ^-Sa ■ vax a v 3 a 



17 

y 



,17 

a A • S A x A ~~ A ' x A x A a A 



,3 ri 43 11 N 4 ( ,. 

^2 A + 28^ ' VA + 14°^ ' 35 A) 2 ^^ 



n STF 



In addition to the spin-orbit multipoles, we need the lowest order contributions to l % i k and J 13 from 
point-mass energy-momentum tensor. These multipoles are given by 



rijk X ^ {^jfe) 

(PM) ~ Z^ mAX A ' 



J ( p M) = ^2m A {x A x v A ) {l x 



A 



(4.8) 
the usual 

(4.9) 
(4.10) 



and contribute to the wave forms at 0.5PN order. 



B. Two-body case 



We now specialize to the case of two bodies of mass mi and vri2 in an arbitrary (possibly precessing) orbit, and 
express the multipoles in terms of the relative coordinate x (whose origin is at body 2). 
It is convenient to use the mass parameters 

M 



mi + vri2, (4.11a) 
mim 2 /M 2 , (4.11b) 
(mi-m 2 )/M. (4.11c) 

It is also convenient to use the dimensionless, symmetrized spin parameters introduced by Will and Wiseman fet, 



V 
A 



Xa 2 



1 /Si 



Xs ~ 2 



1 /Si 



S2 

™2 



We eliminate the potentials and accelerations in (4.; 



"1 "'2. 

with the well-known Newtonian expressions 



OA = 



—M 



-x A , 



(4.12a) 
(4.12b) 

(4.13) 
(4.14) 



4 Although we do not consider spin precession in this paper, the JV-body and two-body multipoles we present are (instanta- 
neously) valid even in the precessing case. One simply has to put in the spin vectors as (slowly varying) functions of time — which 
is easier said than done. 
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where B ^ A and r — \x\. 

The spins of the bodies introduce a correction to the relation between xi, ^2, and the relative coordinate x, 



Xi = x 



X 2 = X 



a ( 2 M 

nA ( v 2 

M 2 V r 



mi 1 



MlJO X (Xa + Axs), 
-Mtjvx (Xa + Ax s ). 



(4.15a) 
(4.15b) 



(Compare Eq. (3.13) of Ref. E3] and Eq. (Fll) of Ref. II , where the missing factor of r\ in the latter is a typographical 
error.) This correction is 1.5 PN order; therefore it enters our 2PN calculation through contributions from the 0.5PN 
(point-mass) multipoles (|4.9|). 



Applying the tra nsfo r matio n (4^158) and including the contributions from the point-mass multipoles, we find the 
two-body forms of ([hj), (pL6j) , and (H-8|) to be 



J^= -MW(Xa-A Xs )^ k *\ 



(4.16) 



f 3 1 STF 

pi k = -M V Ax^ + M 2 v 2 l -[« x ( Xa - 5A Xs )]V fe - 3[a: x ( Xo - A Xs )]V fe \ , 



(4.17) 



3 3 / M\ 

J« = - M A(x x »)<V> + -M 2 v ( Xa + A Xs )<V> + -M 2 tjA U 2 - —J (A Xa + x.)<V> 



~mV 

28 ' 



23y - 13t> 2 ) Xa + A (47 y - 141 ) \ , 
M 



+2(0; • w)(15xa + 13A X .)V + 2-=-[x • (29 Xa - A Xs )]x y - 4[x • (5 Xa + 9A Xs )K 



-4[« • (3 Xa - 3lA Xs )]xV 



STF 



(4.18) 



The first terms in (4.17) and ( 4.1 8| ) are the lowest order non-spin terms. The second term in (4.18) is easily verified 
as identical to the lowest order (1PN) spin-orbit contribution obtained by Kidder [Eq. (3.20a) of Ref. []l2| l. 



C. Quasicircular orbits 



We now specialize to the case where the two bodies orbit each other in a circular trajectory, which adiabatically 
shrinks (inspirals) under radiation reaction. For spinning bodies, this is only possible when both spin vectors are 
parallel or antiparallel to the orbital angular momentum, eliminating spin-orbit precession. 

We express our results in terms of the or thonorma l vectors n = x/r, A = v/v, and z = nx \ (parallel to all angular 
momenta). The majority of the terms in (4.16-4.18) vanish, and we are left with the greatly simplified expressions 



r^ = ^M 2 V 2 r 3 (xa~A Xs )n^ k /K 



(4.19) 



I» k = -M V Ar 3 n^ + -M 2 rfvr 2 [( Xa - 5A Xs )n« fe + 2( Xa - A Xs W^ k 



STF 



(4.20) 



J rj = -Mr)Ar 2 vn {i z j) + |mV(x + A Xs )n ( V> + ^-M 3 r) 2 (5x a ~ 47 A Xs )n {t z 3) . 



(4.21) 
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V. WAVE FORM 



Having evaluated the necessary radiative multipoles, we obtain the gravitational wave form 



h %1 



R^ 



e=2 



_1 pjL—2 ]\fL-2 



c Pl(i J3)pL-2 jyqL- 



(5.1) 



where TT denotes the transverse traceless projection, (^) denotes the £th time derivative, and N is the unit vector 
pointing toward the observer [cf. Eq. (E5a) of Ref. Q]. 

We evaluate the time derivatives using the equation of motion for circular orbits, a — —uj 2 x, where the standard 
form of the post-Newtonian expansion is given by 



uj 2 = 



AI 



1 - (3-77)— - 2[A Xa + (1 + v)Xs}— + 0(e 2 ) 



(5.2) 



[cf. Eqs. (7.1) and (F20) of Ref. [§]. The 1.5PN spin-orbit term in the equation of motion [the third term in (|J)] 
contributes to the 2PN sp in-or bit term in h 1 ^ through time derivatives of the 0.5PN point-mass multipoles [cf. the 
first terms in ( l4~2(il ) and ( l4~2ll)1 If we were to consider spin-o rbit precession, time derivatives of spin expressions 
appearing in the multipoles ( 4.16| )-( fI.18|) would al so f actor into (5.1). 

Evaluating the time derivatives in ( |5.l|) , using (5.2) and the identity v — uur to write v and M/r in terms of Mu, 
and collecting all terms of order (Mw) 2 , we obtain 



R 



(A 2 X a + A Xs )(( 



14 i 



(X- N)- —n^X j) {n- N) 

O 

TT 



2 X „, --A 2 XQ -|A Xs )(nx JV) ( V> 



2M?? 2 r N2 
+ _ R ( } 



19x 



a + 37A Xs ) (i(iV ■ A)n« + 1 -{N ■ n)n^) 



-4( Xa + A Xs )(N ■ A)A y + ( X a - A Xs )<{ (n x iV)M 



on 7 

-_(JV • A)(2V • z)(n x iV)^ + -(TV • n)(7V • z)(n x 7V) ( V) 
3 3 

- — (N ■ n)(N -z)(Xx N'fX^ — — (N ■ n)(N ■ X)(z x N)^X j) 
3 3 



TT 



7 



(N ■ nf 



(5.3) 



Finally, we project out the amplitudes of the h + and h x polarizations. This is done by taking 



h x = -(Pi^ +q l p j )h %: > . 



(5.4a) 
(5.4b) 



Note that the TT projection in (5.1) is subsumed in this operation (cf. Sec. VIIA of U). The standard convention |9J 
is to use the unit triad composed of N (the direction to the observer), p (pointing from the descending node of the 
orbit to the ascending node), and q = N x p. The orbital phase ip is measured from the ascending node, and the 
orbital inclination angle i is given by cosi = z ■ N. Thus we have 



n = p cos ip + (q cos i + N sin i) sin ip, 
X = p sin ip + (q cos i + N sin i) cos ip, 
z = q sin i + N cos i. 



(5.5a) 
(5.5b) 
(5.5c) 



Following |p2| , we organize the amplitude contributions of the wave- form polarizations according to post-Newtonian 
order and physical origin as 



2Af?7 
R 



4% + ■ ■ ■ + x 3 / 2 </ x 2) + x^ 2 H^ SO) + x 2 H ( ^ x + x 2 ij| 2 ; x SO) + x 2 H^ S) 



(5.6) 







where x = (Mw) 2 / 3 . The 2PN spin-orbit contributions to the wave-form polarizations are given by 



H 



(2, SO) 



{ [(109 + 15c 2 ) X a + 7 (1 + 3c 2 ) A Xs ] V + \ (1 + <?) (.Xa + A Xs ) } sin i cos i/j 
+ {^[(H + 5c 2 )xa + (l + 7c 2 ) A Xs ]v- ^(l + c 2 )( Xa + Ax s )}sini cos3?/> 



(5.7a) 



H 



(2, SO) 



{ [(127 - 3c 2 ) X a + (25 + 3c 2 ) A Xs ]v + \{Xa + A Xs )}c sin* sin^ 
+ {~ [(19 - 3c 2 ) X a + (5 + 3c 2 ) A Xs ]v~ ~(Xa + A Xs )}c sini sin3V>, 



(5.7b) 



where c = cosi. 

It is clear that these 2PN contributions to the wave-form amplitudes do not contribute to the radiation reaction 
at 2PN order because the harmonics (u> and 3w) average to zero when beat against the "Newtonian" terms -H+ x at 
frequency 2lo [cf. Eqs. (3,4) of Ref. §§]]. 



VI. SUMMARY 



We have calculated all the non-precessional 2PN spin-orbit effects on the gravitational wave forms of compact 
bodies in quasicircular orbit, and have shown that there is no spin-orbit radiation reaction effect at 2PN order. Our 
calculation was greatly simplified over previous spinning-body post-Newtonian efforts |l3| , |l2| ] by the use of a ^-function 
energy-momentum tensor for spinning particles. We have presented the wave-form polarizations in "ready-to-use" 
form (cf. |22|). 



Note that terms of 0(rf) contribute significantly to (5.7). These are the terms that could not have been obtained 
by the black-hole perturbation approach. Their presence leads us to expect that the 2.5PN radiation reaction will 
also contain significant terms of 0(r/) which are not found in ]l7| , p3| , po| . 

In this paper, we treated the bodies only to linear order in their spins (i.e., considered only spin-orbit effects). Spin- 
spin effects can be treated by a more complicated calculation using the MST energy-momentum tensor. However, 
because spin-spin effects appear at 2PN order, consistency would require that one also include the effects of the bodies' 
quadrupole moments ]24|. 
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APPENDIX: 

In this appendix, we derive the 1.5PN spin-orbit multipoles I %q^ and obtained in p2| 



(SO) 



The general expressions are quite simple, requiring (4.1) and (4.2) only to lowest order, 

jiJ _ / j3 ij ( T 00 ,rpa 

J (so) - a xx ^(so) + J (: 



(SO) 



4 d [ j3 ijkarpQa 
jijk _ ab(i / ,3 jk)arp0b 

(so) — ' ' ■ << 



(SO)- 



(Al) 
(A2) 



Using (2.8) and (3.5) only to lowest order, it is straightforward to evaluate these expressions in the iV-body case as 
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I lio) = J2[ 4x A( v Axs A y 



il{4(^x^}] STP , 



J, 



ijfe _ 



(SO) 



_« cfeiSTF 



(A3) 
(A4) 



Using the transformation (4.15a) to the relative coordinate, we find the two-body forms of the multipoles to be 



jg Q) = |A/V[2z> x X s) j - v\x x X sV] STF , 



jf 0) =AMWw 3 x k s r F - 



(A5) 
(A6) 



These results agree with those obtained using the fluid body energy-momentum tensor |lz 
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